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INTRODUCTION 


Work  on  optima  space-time  filtering  -at-UBL^has  consisted  of  an 
analog  computer  simulation  of  a narrow-band  Mermoz  processor  which  was 
optimum  only  at.  its  center  frequency,  and  of  single  frequency  array 
gain  studies.^  It  is  now  planned  to  simulate  the  Mermoz  processor 
in  non-real  time  on  a general  purpose  computer  using  digital  filters 
and  data  from  a bottom-mounted,  vertical  array.  This  simulation  will 
differ  from  earlier  work  in  three  ways:  (a)  real  sea  noise  will  be 
used  instead  of  electronic  noise  generators,  (b)  the  system  will  be 
optimum  at  more  than  one  frequency  in  the  band,  and  (c)  the  system 
will  be  digital  instead  of  analog.  The  purpose  of  this  report  is  to 
present  the  results  of  a theoretical  study  undertaken  to  predict  the 
perfo'-aance  of  the  simulated  optimum  processor. 


Oefinitlons 


Consider  the  system  shown  in  Figure  1,  where  the  ^2' 3 are 
linear  filters,  s(f)  » F £s(t)]  , F denotes  Fourier  transform,  * 

dc-otes  complex  conjugate,  and  T.  = time  for  a plane  wave  to 
1 from  phone  #1  to  phone  # i for  a given  direction  of  look. 
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Figure  1 

The  conventional  tine- shift -end-sun  processor  results  when  the 
are  given  by  the  following  matrix  equation  at  each  frequency 


V = 


where: 


A'-  IT*  T ,t 

I ( 
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The  optimum  processor  i s obtained  when  the  [ iK  J 

Z =Q"V 


are  given  by: 


where: 


0 = 


QJf) 


0..(- 1) 


0i>(fJ=F  |e  k(*)H,(t*r ,]}  - cro«.-.p.ctr«l 


E [ ] 


denotes  expected  value 


A slope  of  -6  db/ octave  will  be  assumed  for  the  noise  sprectrum 
at  each  receiver. 


QA*  J = ( "f'°-  J for  a11  1 


Since  the  actual  data  will  come  from  a fixed  array,  a surface  generated 
ambient  noise  model  will  be  used  to  specify  the  cross-spectral 
densities.  For  vertical  receivers  the  result  ia: 
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r, 

l~~\  “>  ' / 5 I *J 

'va  .»  = t-  U 1 


LV  i-..  >*  / ?U^7  — / J 


o - 0 

v->  4»  MU. 


where: 


2 Trf  Jlt 


ci1K  * distance  between  1th  and  J*h 
receiver* 


c = velocity  of  sound 


For  the  purpose  of  simplification,  the  simulation  will  be  done 
with  only  two  receivers  and  broadside  steering.  A CW  pulse  will  be 
used  for  the  signal,  i.e., 


S (*)  = ^ 


I A COJ  c il  tc  f , 0 ■<  f « T, 

C*  /O’flirrvvijt’ 


The  top  receiver  will  be  labeled  #1. 
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The  criterion  for  comparing  the  optimum  processor  to  a conventional 
processor  will  be  output  signal-to-noise  ratio.  If  y(t)  ■=  (t) 

for  ar.  Input  of  signal  alone  (see  Figure  l),  and  y (t)  - r (t) 
for  an  input  of  noise  alone , then  the  output  signal-to-noise  ratio 
is  defined  as: 


Equ.  3 


er 


L LJL 


Discussion 

The  amplitude  spectrum  of  the  signal  is: 

e All  e J A'l*  r (Li 

0 ' ' L 1 /J  2 rr  (f  ♦ i. 


AL 

2 


7r(f-f.)T\ 

e ' Tr  q - 

Tr  (+  - fc  JT, 


Essentially  all  of  the  signal  energy  is  contained  in  the  band  defined 
by  the  *ero  crossings  of  S(f)  which  are  nearest  f0  , i.e.,  the 
points  fc  - 4-  and  f0  +- . Hence  there  is  no  need  to  design 
^e  { Z,}  filters  outside  'this  band.  In  other  words,  the  {£,] 
will  be  band-pass  filters  with  arttaetlc  swan  frequency  f*  and 
b'vndv’dth  Since  the  optiimmi  filters  are  band-limited, 

their  impulse  responses  extend  from  - <■—  to  + —>  . Obviously, 

If  the  inpulse  responses  are  to  be  convolved  with  a finite  amount  of 
'em!  data,  they  must  be  truncated.  This  will  be  accomplished  by 
multiplying  each  impulse  response  by  a Burning  window 
flrivr  by: 


'PET. 


USL  Tech.  Memo. 
2242-217-67 


At  this  point  it  is  necessary  to  specify  two  constants  which  will 
be  used  when  real  data  is  recorded.  They  are: 

fo  IS  0 Hi 

T*  M . 1 i > i > 

H 


Thus , the  i_  '£  . } oust  be  detemined  over  the  band  675  Hz  to 
825  Hz.  A resolution  of  10  Hz  will  be  used.  The  following  figure 
shows  the  amplitude  response  of  the  opt  ins  Z.  along  with  the 
approximation  to  it  corresponding  to  T - 4o  as  for  the  Impulse 
response. 
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mu  IN  H7 

3 


The  optimum  and  approximate  phase  functions  are  as  follows 


■r  VI  .(..‘■Vi 


Figure  4 

Clearly  the  approx latt lone  are  good.  It  should  be  Mentioned  that 


The  impulse  response  corresponding  to  the  approximations  of 
Figure  3 and  4 is  shown  In  Figure  5.  The  amplitude  has  been  noneallted 

~o  unity. 


P'** 
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i in  secsixio1 1 


Figure  *> 

Saaplea  of  t , (t ) will  be  used  as  filter  weights  in  the  convolution 
process,  and  it  has  been  determined  that  an  8 ICC  sampling  rate  should 
be  sufficiently  high  to  accurately  represent  r , (t) . Thus  there  will  be 
a total  nuhber  of  filter  weights  for  c,(t)  equal  to 
(0000  samples/sec )(80  ms)  ♦ 1 - 641. 
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Derivation  of  Input  end  Output  81gnal-to-lolae  Ratios 

The  Input  signal -to- noise  ratio  will  be  defined  act 


(s/l ) In  * average  signal  power  at  one  receiver 
average  nolee  power  at  one  receiver 


i 

T 


(W|N  - 


df 


I o- 


;df 


Outfit  elgnal-to-nolse  ratio  has  bean  defined  by  Sqn.  3,  and  vlU  be 
repeated  for  convenient  reference. 


(*/h) 

°-'t  (*; 


An  expression  for 


cr-(f;-  (<ji 


- s(i) 


can  be  derived  in  the  following  Manner. 

_ . 2 n f Tj 


.if  -i  r i..7TT* 

I <r(+J  - «r(fj  eJ  df 


H,+  ) 
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5 ( m 


+ 2.J  df 


An  expression  for  E.\f't(^)\  is  derived  by  starting  with  the  power 
spectna  of  the  noise  after  the  suaser,  which  Is: 


Nr  = V } i.  z. 


t?., 
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Afte-  the  Batched  filter  the  noise  power  spec  truss  is: 


Therefore  the  variance  of  the  output  is: 


df 
J . 


Substitution  of  Squ's  6 and  7 into  Equ.  3 gives: 


Equ.  8 


(5/\Lf 


-5  ( f ) J [i? , + 2 ..  d 4 


N£j(f.)  d ■f 


Results : As  mentioned  earlier,  the  criterion  that  was  chosen 

for  comparing  the  optimum  and  conventional  processors  was  output 
s .gnal-to-noise  ratio.  Therefore  a computer  program  was  written  to 
evaluate  Equ.  8.  It  should  be  mentioned  that  Squ.  8 is  perfectly 
general,  and  therefore  is  valid  for  both  optimum  and  conventional 
(Z4  } . All  parameters  were  chor.en  so  as  to  accurately  model 

actual  conditions  under  which  data  will  be  recorded.  These  conditions 
are: 

(s/N)  in  - -20db 

T s * 13-3  »s 

d = receiver  separation  * 2 ft. 

M = number  of  receivers  * 2 
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Direction  of  look  = broadside 
Noise  spectral  slope  * -6  db/ octave 
Bandwidth  * 150  Hr 

Juration  of  opt  i hub  impulse  responses  * 2 T • 60s a. 

Sampling  rate  = 8 KHz 

The  resulting  output  s ignal -to-noise  ratios  were  found  to  be: 

(S/H)  out  * -12.5  db 

conventional 


(8/N)  out  *=  -5.3  db 

optima 

Thus  overall  conventional  processing  gain  is  7.5  db,  as  compared  to 
14.7  db  for  the  optima  systea.  The  difference  of  about  7 db  between 
these  two  figures  Is  attributable  primarily  to  the  optimum  {Z.,}-  , 

as  one  discovers  if  he  computes  and  casernes  array  gains  for  the 
frequencies  and  spacing  of  receivers  in  this  systea. 

It  is  not  expected  that  processing  gains  of  7.5  db  and  14.7  db 
will  be  obtained  when  real  data  is  filtered  in  the  future  simulation 
because  the  preceding  calculations  have  been  made  under  the  assumption 
of  signal  known  exactly,  whereas  in  practice  the  amplitude  and  time 
of  arrival  will  not  be  exactly  known,  and  there  may  be  slight  frequency 
spreading  of  the  pulse  due  to  a required  reflection  from  the  surface. 
However,  the  difference  between  optimum  and  conventional  processing 
should  be  roughly  the  same  as  calculated  above  since  the  Inexact 
knev'.edge  of  the  s gnal  should  degrade  both  processors  to  about  the 
sar.**  extent. 


12 


USL  Tech.  Memo. 
2242-217-67 


REFERENCES 


1.  Kennedy,  Robert  M.  and  Robert  O’Neill,  "Adaptive  Beamforming  for 
Sonar  Signal  Processing",  USL  Technical  Memorandum  No.  914-184-66, 

June  29,  1966. 

2*  Becker,  Clair  J.  and  Benjamin  F.  Cron,  "Optimum  Array  Oaln  for 
Di-ectional  Noise",  USL  Report  No.  656,  Oct.  5,  1965- 

3-  Cron,  Benjamin  F. , Norbert  P.  Fisch,  and  Robert  L.  Shaffer,  "Array 
Cain  for  the  Case  of  Horizontal  Noise",  USL  Report  No.  731,  April 

1966. 

4.  Cron,  Benjamin  F.  and  Ronald  R.  Kneipfer,  "The  Effect  of  Spatially 
'.’ncorrelated  Noise  on  the  Maximum  Array  Gain  Processor",  USL 
'“-'chnical  Memorandum  No.  2242-324-66,  Dec.  15,  1966. 

5.  Cron,  Benjamin  F.,  Beverly  C.  Hassell,  and  Frank  J.  Keltonic, 

Addendum,  Spatial -Correlation  Functions  for  Various  Noise  Models", 
J.  Acoust.  Soc.  Am.,  Vol.  38,  p.  885,  Nov.  1965. 


13 


